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It was reported earlier [Phys. Rev. Lett. 106, 056401 (2011)] that the skutterudite structure 
compound CoSbs displays a unique band structure with a topological transition versus a symmetry- 
preserving sublattice (Sb) displacement very near the structural ground state. The transition is 
through a massless Dirac-Weyl semimetal, point Fermi surface phase which is unique in that (1) it 
appears in a three dimensional crystal, (2) the band critical point occurs at fc=0, and (3) linear bands 
are degenerate with conventional (massive) bands at the critical point (before inclusion of spin-orbit 
coupling). Further interest arises because the critical point separates a conventional (trivial) phase 
from a topological phase. In the native cubic structure this is a zero-gap topological semimetal; we 
show how spin-orbit coupling and uniaxial strain converts the system to a topological insulator (TI) . 
We also analyze the origin of the linear band in this class of materials, which is the characteristic that 
makes them potentially useful in thermoelectric applications or possibly as transparent conductors. 
We characterize the formal charge as Co^ d*, consistent with the gap, with its 3 site symmetry, 
and with its lack of moment. The Sb states are characterized as px (separately, py) cr-bonded S&4 
ring states occupied and the corresponding antibonding states empty. The remaining (locally) pz 
orbitals form molecular orbitals with definite parity centered on the empty 2a site in the skutterudite 
structure. Eight such orbitals must be occupied; the one giving the linear band is an odd orbital 
singlet A2u at the zone center. We observe that the provocative linearity of the band within the 
gap is a consequence of the aforementioned near-degeneracy, which is also responsible for the small 
band gap. 

PACS numbers: 



I. MOTIVATION AND INTRODUCTION 

Recent years have seen an explosion of interest in 
two areas that, while distinct in themselves, have 
also found some areas of overlap. One is the situa- 
tion where Dirac-Weyl linear bands emanate from a 
point Fermi surface (FS), the so-called Dirac point. 
While linear bands around symmetry points are not 
uncommon in crystal structures, it is uncommon 
that such points can determine the Fermi energy 
(Ei?), which is what happens in the celebrated case 
of graphene. It is necessary that there is a gap 
throughout the Brillouin zone except for the touch- 
ing bands. Point FSs occur also in the case of con- 
ventional (massive) zero-gap semiconductors ji and 
give rise to properties that have been studied in some 
detail. The case of the Dirac-Weyl semimetal (point 
Fermi surface) systems is quite different j^i^ leading 
to a great deal of new phenomenology. 

The other new and active area relevant to this pa- 
per is that of topological insulators (TIs), in which 
the Brillouin zone integral of a certain 'gauge field' 
derived from the fc-dependence of the periodic part 
of the Bloch states of the crystal is non-vanishing^Ti^ 
These integer- valued topological invariants delineate 



TIs from conventional ('trivial') insulators, with dis- 
tinctive properties that are manifested primarily in 
edge states i^ii"— Identification of TIs, and study of 
their properties and the critical point that separates 
the phases, is highlighting new basic physics and sev- 
eral potential applications. 

Recently we reportedi^ on the skutterudite struc- 
ture compound CoSba, which combines and in- 
terrelates several of these properties (Dirac-Weyl 
semimetal at fc=0; critical point of degeneracy with 
massive bands; transition to topological insulator) 
in a unique way. The skutterudite structure, which 
is illustrated in Fig. [T]and is discussed in Sec. II, is 
critical to the electronic behavior in CoSb3 though it 
may not be necessary for the unique type of critical 
point that arises. The relevant portion of the band 
structure, in a very small volume of the Brillouin 
zone (BZ) centered at P, is unusually simple com- 
pared to most of the TIs that have been discovered, 
viz. the Bi2Se3 classyii, the HgTe clasSfi^ and the 
three dimensional system HgCr2Se4J^ 

The electronic structure of skutterudites was stud- 
ied early on by Singh and Pickett who focused on 
the peculiar valence band which was linear (except 
exceedingly near P) and whose linearity extended 
surprisingly far out into the BZ. It became clear 



that the Unearity of the band was responsible for 
the large thermopower— that is potentially useful 
in thermoelectric applications J^iii The origin of the 
linearity was however not identified, but it leads to 
peculiar consequences the density of states varies 
as £^ near the band edge rather than the usual 
three dimensional (3D) form y^; as a consequence 
the carrier density scales differently with Fermi en- 
ergy ep; the inverse mass tensor VV^fe is entirely 
off-diagonal corresponding to an "infinite" transport 
mass; the cyclotron mass differs from conventional 
3D behavior, etc. As the Sb sublattice is shifted in a 
symmetry-preserving way^ the small semiconduct- 
ing gap of CoSba closes at a critical point, giving 
rise to an unusual occurrence: a Dirac-Weyl point 
in a 3D solid at the F point. This transition point 
sets the stage for a conventional insulator to topo- 
logical insulator transition, although the specific be- 
havior of the system at this point was not spelled 
out in our earlier work. The only other topological 
insulator phases in skutterudite materials that we 
are aware of are the examples by Jung et alJ^ in 
CeOs4Pni2 filled skutterudites. 

This paper is organized as follows. In Section 
II the skutterudite structure and its relation to the 
symmetric perovskite (ReOa) structure is reviewed. 
The computational methods are outlined in Section 
III. Section IV discusses the evolution of the band 
structure during a perovskite-to-skutterudite trans- 
formation, including the development of the gap and 
the appearance of the linear band. Section V is de- 
voted to constructing a microscopic but transparent 
understanding of the electronic structure, including 
the development of the energy gap and the character 
of the peculiar linear band and the lower conduction 
band triplet that it interacts with near, at, and be- 
yond the critical point. In Section VI we provide 
details of the band inversion leading to the topo- 
logical transition, and analyze the anisotropy at the 
critical point. A brief summary is provided in Sec. 
VII. 



II. STRUCTURE AND RELATION TO 
PEROVSKITE 

It is useful for the purposes of this paper to con- 
sider the skutterudite structure to be a strongly dis- 
torted perovskite ATPua structure in which the A 
atom is missing {i.e. the ReOa structure) and the 
interconnected TPng octahedra arc rotated substan- 
tially (T — transition metal; Pn — pnictogcn). The 
skutterudite structurer^"— pictured in Fig. [U has 
space group Im3 (#204) and a body-centered cubic 
(bcc) Bravais lattice, and is comprised of a bcc rep- 




FIG. 1: (Color online) Crystal structure of skutterudite 
minerals {viz. CoSbs) with space group ImS (#204), 
which includes inversion. The Co site (small pink [dark 
grey] sphere) is octahedrally coordinated to Sb atoms 
(small yellow [light gray] spheres). Each Sb atom con- 
nects two octahedra, as in the perovskite structure which 
has the same connectivity of octahedra. The large (blue) 
sphere denotes a large open site that is unoccupied in 
CoSbs but is occupied in filled skutterudites (see for ex- 
ample Ref. [3l). The geometrical solid (center of figure) 
provides an indication of the volume and shape of the 
empty region. 



etition of four formula units (f.u.) when expressed 
as TPn^. The pnictide {Pn) atoms form bonded 
units (nearly square but commonly designated as 
"rings" ) which are not required by local environment 
or overall symmetry to be truly square. The three 
Pn4 squares in the primitive cell are oriented per- 
pendicular to the coordinate axes. Transition metal 
(r) atoms (usually in the Co column) lie in six of 
the subcubes of the large cube of lattice constant a; 
the other two subcubes (octants) are empty. The 
structure has inversion symmetry and is symmor- 
phic, with 24 point group operations. The cubic op- 
eration that is missing is reflection in (110) planes. 
The related filled skutterudites XT4Pni2 have an 
atom X incorporated into the large 2a site of 3m 
symmetry^^ that remains empty in the compounds 
that we discuss. 

Figure [5] illustrates two aspects of the structure. 
The top panel reveals that, although the Sbe octa- 
hedra are rotated considerably (and strained) com- 
pared to the perovskite structure, a great deal of reg- 
ularity is retained. The bottom panel gives a picture 
of the relationship between the Co atoms and the 
six Sb4 rings that it is coordinated with (only three 
mutually orthogonal rings are pictured). A natural 
local coordinate system for describing an Sb4 unit is 
to let it lie in the x — y plane; then the Pz orbitals 
we will discuss will be oriented perpendicular to the 
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FIG. 2: (Color online) Top: view down a cubic axis of 
CoSba, illustrating the very regular appearing positions 
of the Sb4 rings, drawn as connected by bonds. Note 
that one of these rings looks like a dumbbell when viewed 
edge on. The (dark) Co atom is pictured without bonds. 
Bottom: an indication of the environment of the Co atom 
(smaller sphere) ; only three of the six attached Sb4 rings 
are pictured. In this structure, all Sb sites are equivalent 
as are all Co sites. 



plane of the Sb unit. The Px,Py orbitals are nearly 
symmetry related and naturally lead to bonding and 
antibonding molecular orbitals, and hence are very 
similar in character, while the nonbonding pz or- 
bitals have distinctive character. 

If the origin is chosen so the large empty sites are 
centered at and (|, |, |), then the rings are 

centered in each of the other subcubes (3:, |, |) etc. 
A ring oriented perpendicular to z is neighbored by 
rings above and below (±z directions) oriented per- 
pendicular to (say) y, by rings in the ±i directions 
perpendicular to x, and neighbored by the empty 
sites along the ±y directions. A symmetric com- 
bination of the four pz orbitals on a ring, call it 
P^, is orthogonal to the orbitals on neighboring 
rings, so if dispersion is governed by inter-ring hop- 
ping (rather than through Co atoms) some rather 
flat Sb pz bands should result. Such behavior is seen 
in fatbands plots (sec below). 

As mentioned, the skutterudite structure is re- 



lated to the perovskite structure DTPus (□ denotes 
an empty A site). Beginning from perovskite, a ro- 
tation of the octahcdra keeping the Pn atoms along 
the cube faces results in the formation of the (nearly 
square) Pn4 rings, and the Pn octahcdra become dis- 
torted and less identifiable as a structural feature. 
The transformation is, in terms of the internal coor- 
dinates u and w, 

u'{s) = -^+s{u-]); v'{s) = -^+s{v-^). (1) 

The transformation path, from perovskite for s=0 to 
the observed structure for s=l, is pictured in Fig. 1 
of Ref. [l^. Below we make use of this transformation 
to follow the opening of the (pseudo)gap between 
occupied and unoccupied states. 



III. COMPUTATIONAL METHODS 

Two all-electron full-potential codes, FPLO-92^ 
and WIEN2k^ based on the augmented plane 
wave-|-local orbitals (APW-l-lo) method;^ have been 
used in these calculations, with consistent results. 
The Brillouin zone was sampled with regular 12 x 
12x12 fc-mesh during self-consistency. For WIEN2k, 
the basis size was determined by ^mt^max= 7. 
Atomic radii used were 2.50 a.u. for Co and 2.23 a.u. 
for Sb, and the Perdew-Wang forrr^ of exchange- 
correlation functional was used. The experimental 
lattice parameters and atomic positions were taken 
from Ref. [IMi 



IV. PEROVSKITE TO SKUTTERUDITE 
TRANSFORMATION 

The development of the electronic spectrum as the 
crystal is distorted from perovskite to the skutteru- 
dite structure is pictured in the three panels of Fig. 
[21 For the ideal perovskite the result is a highly 
metallic state, with no gap or pseudogap near the 
Fermi level. The important features only arise near 
the end of the distortion path. For s=0.75 (top panel 
of Fig. [3]) the valence bands have just become dis- 
joint from the conduction band and gap formation is 
imminent. For s=0.90 the gap is well formed and the 
minimum direct gap is emerging at the F point. The 
unusual high velocity band arising from the F point, 
already striking at s=0.75, remains unchanged at 
this point. The valence bands are a mass of indeci- 
pherable spaghetti. 

By s=0.95 (bottom panel) a similarly high ve- 
locity band is emerging from the valence bands to 
the maximum at the F point. At s=1.0 it becomes 



3 



clear that it is a partner of the high velocity conduc- 
tion band, as discussed in our earlier paper— and 
to which we return in the next section. No doubt 
the formation of the gap, i.e. formation of occupied 
bonding bands and unoccupied antibonding bands, 
is behind the stability of the skutterudite structure 
in the Co pnictides class of materials. The fact that 
the gap only opens as the Sb Px and Py bonding- 
antibonding splitting becomes strong supports the 
picture that this tpp^ interaction takes precedence 
over the Co-Sb interactions within the CoSbg unit; 
however, these latter interactions {i.e. the presence 
of the Co atom) are necessary for the gap formation. 
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ELECTRONIC STRUCTURE 



A. The Zintl Viewpoint 

Although in the 'parent' perovskite system the 
structure is that of vertex-connected CoSbe octahe- 
dra, in the strongly distorted skutterudite structure 
the Sb4 rings form a basic structural motif. These 
compounds, viz. CoSba, are sometimes character- 
ized as Zintl materials in which the Sb4 unit balances 
the charge of the Co unit. Within the Zintl picture, 
the charge of the Co™"*" ion must be balanced by a 
(Sb4)(''™/3)- ^nit: there are only 3/4 as many Sb4 
rings as Co atoms. Thus there is tension between the 
Zintl picture, which does not lead to integral formal 
charges on both of the primary units, and the pres- 
ence of the gap that specifies an integral number of 
occupied bonding states. Several paper a^^i^^'^^i^^" — 
have presented results and some analysis for empty 
and for a few filled skutterudites. We have benefited 
from the previous work as we proceed on a deeper 
analysis. 

We first follow the commonly held line of reason- 
ing that the short Sb-Sb distance in the Sb4 ring 
is the most fundamental aspect of the electronic 
structurei^ The Px — Px bonding along an x-axis 
is characterized by a hopping amplitude tpp^ ~ 3 
eV, so the corresponding bonding and antibonding 
states lie at Sp ± tppg-. The separation of Sb atoms 
and the a bonding and antibonding between py or- 
bitals along the y-direction is indistinguishable from 
that of the Px cr case. Indeed, the corresponding 
Per bonding and antibonding character (not shown 
in figures) is centered roughly ±3-4 eV below and 
above the gap. The on-site energies Sp^ ~ Sp^ = Sp 
therefore lie in the vicinity of the gap. This picture 
leads to half- filled Px and py orbitals. 

The Pz orbitals within an Sb4 unit are orthogonal 
to the Px,Py orbitals and couple by a tpp-jr hopping 
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FIG. 3: Transformation of the band structure from 
metallic in the perovskite to semiconducting in the skut- 
terudite, following the (linear) transformation described 
in the text. From the top, the panels show s — 0.75, 
0.90, 0.95. The gap only begins to emerge around s— 
0.75, where the conduction bands are already disjoint 
from the valence bands. The linear band near F emerges 
from the valence bands around s=0.95. 
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FIG. 4: (Color online) Band structure of CoSbs along 
high symmetry directions, with fatbands highlighting the 
Sb pz character. The character is substantial in the -3 to 
-6 eV region, in the 3 to 5 eV region, but most notably 
in the linear band extending through the gap, which is 
primarily pz in character. 

amplitude. This coupling leads to doubly degenerate 
levels at ^tpp^^ relative to Sp. The evidence is that 
there is negligible crystal field splitting of the three 
on-site p orbitals. However, inter-unit — Pz inter- 
action will be stronger than this intra-unit coupling, 
and the Pz orbitals also couple to the Co orbitals. 
The Pz character, provided by the fatbands repre- 
sentation in Fig. 21 is in fact distributed throughout 
the -6 cV to -|-5 eV region and is not reproducible 
by a simple model due to other couplings. The Pz 
projected density of states (DOS) indicates the Pz 
orbitals are at least half-filled, perhaps slightly more. 



B. Complications with the Zintl Viewpoint 

It is readily demonstrated that analysis focusing 
primarily on the Sb4 rings has strong limitations. 
We have performed calculations with the Co atom 
removed, DSba. The Sb p band complex is broad (10 
eV) and featureless. There is no hint of a gap or even 
pseudogap, indicating that the px and Py bonding- 
antibonding feature described above is far from a 
dominant feature. Only when Co is re-inserted into 
the structure does the remarkably clean 1 eV gap 
become carved out around the Fermi level, except 
of course for the peculiar linear band near the zone 
center. The gap therefore arises from the Co Zd mix- 
ing (evidently strongly) with Sb bp orbitals, with the 
linear band emerging from the valence bands being 




FAH NLTAP 



FIG. 5: (Color online) Band structure of CoSbs along 
high symmetry directions, with fatbands highlighting the 
Co 3ci character. This plot reflects a narrow mass of d 
bands throughout the -1 to -2 eV range, more around -3 
eV, and ~4 unoccupied bands in the 0-1 eV range. The 
amount of 3ci occupation is discussed in the text. 



of strong Sb character. 

Figure [S] illustrates the Co 3d character. It lies 
mostly below the gap, in the very narrow -1 to - 
2 eV range. A much smaller amount of character, 
roughly something like four bands, lies immediately 
above the gap in the 0-1 eV range. These two parts 
of the Sd projected DOS are neatly and impressively 
split by the 1 eV gap. The "charge state" (or for- 
mal valence, or oxidation state) of the Co atom is a 
question that can be asked, since this is an insulat- 
ing material. If there arc four unoccupied d bands 
(though this cannot be claimed very conclusively due 
to the strong hybridization) then considering there 
are four Co atoms in the cell, the occupation would 
be d^ Co^+. A nonmagnetic ion with this filling is 
natural for the Co 3 site symmetry. As the CoSbg 
octahedron is distorted strongly in progressing from 
the cubic perovskite structure to the skutterudite 
one, the crystal field levels reduce as t2g — >■ + ^2, 
the latter being the doublet, and Cg — ?> ei -t- 62. One 
unoccupied orbital singlet (both spins) gives the d^ 
occupation. The magnetic quantum number m = 
orbital with respect to the local three-fold axis is the 
natural one for the id holes to occupy. 

Charge balance will then leave the rings as 
(Sb4)~''/'^, a very unsatisfactory result for account- 
ing for the gap: each Sb4 ring would not have an 
integer number of (fully) occupied bands (or molec- 
ular orbitals). Within the picture that the Px and Py 
orbitals are half-filled due to the large tppcr bonding- 
antibonding splitting, this Co charge state would 
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suggest Pz orbitals to be 1/6 more than half- filled, 
a peculiar result and one that would only be com- 
patible with a metallic, rather than semiconducting, 
state. 

The integrated DOS for the Co 3d states in fact 
gives a clear (P occupation. However, with near- 
neutral atoms there might be significant Co 4s, 4p 
occupation, which could bring Co nearer to a neu- 
tral d^s^ configuration. Some Co sp character is in 
fact found in the upper valence bands, though it is 
possible that this represents tails of Sb p orbitals. 
The Co d occupation d^ is clear, however. 



C. Sb 6p - Co 3d Mixing 

It was noted above that when Co is removed from 
the structure, there is no gap nor any indication of 
one, nor any candidate for a potential linear band. 
The most direct interaction of the Sb Pz orbitals on 
the ring is between the Pz orbital on the nearest Sb in 
the ring with the rotated xz±yz orbitals on the eight 
Co neighbors. Sbg molecular orbitals (MOs) formed 
by linear combinations of Pz orbitals can be sorted 
according to whether they are even or odd parity 
with respect to the center of inversion at the Co site. 
The even parity MOs mix with the Co 3d orbitals, 
forming bonding and antibonding pairs tending to 
open a gap; the odd parity MOs do not mix with 
the 3d orbitals. As a result of this coupling, four of 
the five 3d orbitals are lowered and occupied, while 
one is raised and unoccupied, giving the 1 eV gap 
and the d^ occupation. 



D. Band Character near the Gap 

The fatbands representation of Sb pz character 
in Fig. m shows the substantial Pz character of the 
linear band in the upper valence bands near the 
zone center. This pz character is in fact the dom- 
inant character, although various other small con- 
tributions from other valence orbitals are mixed in 
by the low Sb site symmetry. The character of the 
3-fold degenerate conduction band minimum state 
is more complex. These bands have primarily Co 
"i2g" character; the quotes arise from the fact that 
the degeneracies of the cubic t2g and Cg subshells 
are broken by the distortion of the octahedra in the 
skutterudite structure. There is a strongly dz2 flat 
band just above (-1-0.2 eV) the gap, and additional 
strong dx2_y2 character in the 0.2-0.6 eV region. 

The band figures we show below demonstrate that 
the linear band couples to only one of the triplet of 
conduction bands at and near F, i.e. couples to only 



one linear combination of the triplet states. It does 
not couple to the Co 3d orbitals. 




FIG. 6: Surface plot of the density arising from the 
linear band in CoSba near fc=0. The closed surface sur- 
rounds the unoccupied site in the skutterudite lattice 
and encloses a region of lower density than outside of 
the surface. 



E. A consistent viewpoint: an empty-site 
molecular orbital 

A generalized viewpoint of the most physical or- 
bitals in CoSba is that of four Co^"*" d^ atoms (as 
discussed above), three Sb4 rings with p^- and Py-a 
bonding states filled and antibonding states unoccu- 
pied, and molecular orbitals formed from the twelve 
Pz orbitals and centered on the empty 2a site. This 
site is an inversion center, so the MOs can be clas- 
sified as parity-even or -odd. (Note that these are 
different MOs than were discussed in Sec. V.C.) 
Twelve MOs can be formed, which can be classified 
also by the tetrahedral 2a site symmetry. Charge 
counting requires that eight of these MOs must be 
occupied (with both spins). 

Wc will not attempt to identify the positions and 
dispersions of these MOs, except to note that linear 
band is one of these MOs. Jung et al^ identified 
this MO as ^211, an orbitally nondegenerate state 
with odd parity at fc=0, and therefore not mixing 
with the Co 3d orbitals at (or near) F. The density 
arising from this orbital is of interest. We had earlier 
investigated the possibility that an electron may re- 
side in the large empty 2a site in the skutterudite 
structure. This type of "electride" configuration, 
in which an interstitial electron without a nucleus 
becomes an anion, is well established in molecular 
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solidsi^i^ However, the density at the 2a site is very 
smaU, ruhng out this possibihty. 

An isosurface plot of the density arising from the 
Uncar band near F is shown in Fig. [SI The closed 
surface in the density encloses a region of low density 
centered on the unoccupied site. Exploring the den- 
sity with isosurfaces at various values of density re- 
veals no local maximum at this large interstitial site, 
thus no electride-like character. Lefebvre-Devos et 
al. have presented^ complementary isosurface den- 
sity plots of this band, with the high density regions 
of two types. One region surrounds the empty site 
(with rather complex shape) consistent with its ori- 
gin as a MO comprised of Pz Sb4 ring states, with 
the tetrahedral symmetry of the 2a site. The second 
region is within each CoSbg octahedron, centered on 
either side of the Co atom along the local 3 axis, ap- 
pearing as a three-bladed propeller that arises from 
the lobes of the pz orbitals that lie closer to the Co 
atom. 

Interaction of the Pz states with Co 3d orbitals 
appears to cause a majority of the density of this 
state to border the large empty site. The result is not 
an electride state but rather a large molecular orbital 
(one of twelve in total) centered on the empty site, 
with one of them (A2«) per primitive cell. We have 
previously reported that this state at F is parity- 
odd^ 

Such an orbital has the same hopping amplitude 
t along each nearest neighbor direction on a bcc lat- 
tice, hence its dispersion looks like that of an s or- 
bital, Ek = 8t cos{kxa/2) cos{kya/2) cos{kza/2), and 
does not lead to a linear band at k=0. This shouldn't 
be surprising, as the true linearity at F is acciden- 
tal, due to a degeneracy that has to be tuned. In 
our previous paper^i^ this tuning was accomplished 
by making two on-site energies in the tight-binding 
model identical. 

The "linear band" aspect of the skuttcrudites may 
have been over-interpreted, because it arose in the 
most common member (the mineral skutterudite 
CoAsa) which was studied most heavily initially. In 
a general skutterudite, for example the filled one 
LaFe4Pi2, the same band rises out of the valence 
band complex but is not unduly "linear." The iso- 
valent members XPn3, with X = Co, Rh, Ir, and 
pnictides Pn = As and Sb, all seem to lie not far 
from the critical degeneracy that extends the band 
linearly to fc=:0; however, we know of nothing in 
general about the skutterudite crystal structure that 
otherwise promotes a linear band. 



VI. NONANALYTICITY, ANISOTROPY, 
AND THE TOPOLOGICAL INSULATOR 
PHASE 

A. Band Inversion: Internal strain, spin-orbit 
coupling, and tetragonal strain 

In previous workiS we demonstrated that it is pos- 
sible, with very small Sb sublattice displacement (see 
Sec. II for the linear distortion, characterized by the 
parameter s), for this band to remain linear all the 
way to fc=0, marking a critical point. In that work, 
the progression of the band crossing at the critical 
point was plotted. Due to the 3-fold degeneracy of 
the conducting band edge, the progression versus s 
was from semiconductor for s < 1.19, to the critical 
point at 5=1.020 with linear valence and conduc- 
tion bands degenerate with two quadratic bands, to 
a zero-gap semiconductor for s > 1.021 due to a 
symmetry-determined degeneracy, even when spin- 
orbit coupling (SOC) is included. Due to the band 
inversion and the character of the states at F, the 
transition is from trivial insulator to one with topo- 
logical character. Because the degeneracy at the 
point Fermi surface (at F) remains, the latter state 
is a topological zero-gap semiconductor rather than 
a topological insulator. 

We now demonstrate that a simple tetragonal 
distortion, such as might result from growth on a 
substrate with small lattice mismatch, does pro- 
duce a topological insulator (TI) but with an un- 
expected placement of the gap. The progression of 
the band structure very near F is shown in Figure 
[7] for strain c/a=1.01. The gap before the critical 
point is reached always has quadratic bands for suf- 
ficiently small k. At the critical point s = 1.020 
(when the strain is imposed) leads, unexpectedly, to 
linear bands only in certain of the high symmetry di- 
rections (such as the < 111 > direction, toward the 
P point). In other high symmetry directions (from 
F toward H is pictured) the bands are quadratic and 
the upper one is quite flat (heavy mass). We return 
to this point in the next subsection. 

Beyond the critical point, the four bands become 
non-degenerate once more; the lower (linear) band 
has crossed the lowest of the three conduction bands 
and acquired mixed character. The coupling be- 
tween bands is clearly much smaller along the axes 
(F-H direction) than along the other symmetry lines, 
and the gap occurs along this direction. Because of 
the heavy mass of the lower conduction bands, the 
gap is quite small, being of the order of 1 meV for 
the tetragonal strain that is pictured in Fig. [T] 
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< 100 > directions can be constructed as 



Piooik) 



Y[\k-k,\, 



(2) 



where {kj} are the six < 100 > unit vectors. The de- 
nominator normahzes the polynomial to unity along 
each of the < 111 > directions; km lies along 
the < 111 > direction (any of them). The anal- 
ogous polynomial Piio{k) that vanishes along each 
< 110 > direction and is normalized to unity along 
each < 111 > direction can be constructed. Then 
the peculiar band pair is given by 



e{k) = ±v\k\PwoCk)Piio{k) + 



2m, 2m„ 



(3) 



where nis , vria serve as symmetric and antisymmetric 
masses. 



Isovalent Skutterudite Antimonides 



FIG. 7: Transformation of the band structure very 
near k=Q from insulating at s= 1.019 (top) through the 
semimetallic critical point at s=1.020 to the inverted 
band insulator at s= 1.022 (bottom). H/50 and P/50 
denotes the point 2% of the way from F to H and P, re- 
spectively. The horizontal line lies in the gap (or at the 
Fermi level) in each panel. In the bottom panel the gap 
has moved out along the F-H line but is non- vanishing. 
The much larger band repulsion along F-P produces the 
unique band dispersion at the critical point. 



B. Anisotropy at the Critical Point 



Inspecting the bands along the three cubic high 
symmetry directions, a second peculiar feature be- 
comes evident. At the critical point, the "linear 
band" disperses linearly along the < 111 > di- 
rections, while they are quadratic (with distinct 
masses) along both the < 100 > and < 110 > di- 
rections. The directional dependence as well as the 
magnitude dependence of the wavevector is anoma- 
lous. The form of the dispersion for small \k\ can 
be constructed. A polynomial in the direction of 
k (k) with cubic symmetry that vanishes along the 



We have studied the band structures of anti- 
monidc sisters of CoSba, based on the transition 
metal atoms Rh and Ir. Band structures, or some- 
times only gap values, of these compounds have been 
reported earlier j^^^i^^ The band structures around 
the gap at fc=0 differ somewhat in the literature, 
due to sensitivity to structure (and that some use 
relaxed coordinates while others use the measured 
structure), to the level of convergence of the calcu- 
lation, and the exchange-correlation functional. 

The difference in lattice constant is compensated 
very closely by the effect of difference in interatomic 
distances on hopping amplitudes, and the resulting 
bands are very similar. The single feature of inter- 
est is the magnitude of the gap (and whether posi- 
tive or negative, the latter giving an inverted band 
structure) and the strength of SOC effects. As men- 
tioned above, band inversion alone does not in it- 
self give rise to a topological insulator state in these 
compounds because even after SOC is included, a 
two-fold degeneracy at F is lowest and the system is 
a point Fermi surface, zero-gap semiconductor. Such 
cases are referred to as topological semimetals. How- 
ever as shown above, the topological semimetal can 
be driven into the TI phase by tetragonal distortion 
that lifts the final degeneracy. 

RhSb-^. For a range of volumes around the ex- 
perimental one (a= 9.23 A), this compound has an 
inverted band structure, hence it is a point Fermi 
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surface, zero-gap topological semimetal. In the in- 
verted band structure before considering SOC, the 
singlet lies 23 meV above the triplet that pins the 
Fermi level. Inclusion of SOC leaves the Fermi en- 
ergy pinned to a doublet. A tetragonal distortion 
raises the degeneracy; a topological insulator is ob- 
tained for c/a > 1 in the same way that happens 
for CoSbs although no sublattice strain is required 
in RhSba. 

Generally, reducing volume (applying pressure) in 
these Sb-based skutterudites leads to a larger gap. 
In the case of RhSba , our calculations indicate a gap 
never actually opens within cubic symmetry, for pe- 
ripheral reasons. When the band inversion at F is 
undone and a gap at F appears at sufficiently high 
pressure (a= 9.06 A), the material is metalhc due 
to other (Rh 4d) bands having become lowered in 
energy, crossing the Fermi level along the F-N direc- 
tion. 

IrSbs- The larger SOC for heavier atoms led us 
to consider this 5d compound. However, the linear 
band and several others have practically no metal 
atom character and are not affected by its spin-orbit 
strength. We obtain, using the same computational 
methods as for the other two compounds, that like 
CoSbs, IrSba has a ^80 meV gap and therefore no 
topological character at the experimental equilib- 
rium structure. Like CoSba, it can be driven to an 
inverted band structure and hence topological char- 
acter by internal and tetragonal strains. 

VII. SUMMARY 

One goal of this work was to provide a simple but 
faithful picture of the electronic structure of CoSba : 
the origin of the gap and the character of the linear 
band being the most basic. The picture is this. Co 
is in a configuration, leaving four electrons in the 



Co4Sbi2 unit cell to go into other bands. The Px 
and py orbitals on the Sb4 ring form strongly bond- 
ing (occupied) and antibonding (unoccupied) states, 
accounting for two of the three 5p electrons of each 
Sb atom, for a total of twelve occupied cr-bonding 
states (of each spin). The twelve 7r-oriented Pz or- 
bitals couple between units as well as with the Co 
3d orbitals, and form "molecular orbitals" centered 
on the vacant 2a sites in the skutterudites lattice. 
Eight of these MOs are occupied by 16 electrons, 
one from each Sb and one from each Co. The linear 
band is best pictured as arising from the coupling 
and dispersion of these MOs. 

We have also justified our earlier statementi^ that 
CoSba is very near a conventional-to-topological 
transition. A small symmetry-preserving internal 
strain, a small applied tetragonal strain, and spin- 
orbit coupling drives CoSba into a topological insu- 
lator phase. This critical point also marks a point 
of accidental (but tunable, by the strains) degener- 
acy that gives rise to emergence of linear Dirac-Weyl 
bands emanating from fc=0, where they are degen- 
erate with massive bands. At this critical point, the 
dispersion is non-analytic at k~0 and anisotropy is 
extreme, with the mass varying from normal to van- 
ishing. In this respect CoSbs at its critical point 
bears several similarities to the semi-Dirac (Dirac- 
Weyl) behavior encountered^"— in ultrathin layers 
of VO2. 
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